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Abstract 

An aecelerator-driven subcritieal reaetor functions weU in incinerating high-level radiotoxic waste (HLW) as 
well as providing energy. China is on his way to establish such a facility to transmutate the annual 1000 
tons of HLW. A neutronic analysis has been performed for a reference core with a special task of burning 
minor actinides ^"^^Np, ^""^^Am, ^^'^Am and ^^"^Cm. Instant operation parameters are determined, including 
neutron energy spectra, thermal power distribution and transmutation validation. Burning-up analysis is 
carried out to further confirm the incineration efficiency. The core parameters optimized in this work will 
be applied to simulate in-core fuel behavior in future research. 
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1. Introduction 



Nuclear Power is one of the major sources of en- 
ergy in the world. The installed capacity of the 
436 nuclear power reactors operating in the world 
is 370 GWe and these provide nearly 18% of the 
worlds electricity [1]. However, along with the en- 
ergy, a large amount of high level radioactive waste 
is discharged from the spent fuel every year. Reduc- 
tion of burden caused by radioactive waste man- 
agement is one of the most critical issues for the 
sustainable utilization of nuclear power. Without 
efficient transmutation, time evolution of the ra- 
diotoxic waste is shown in Figl. The toxicity de- 
pletion time is in the scale of million years [2]. An 
accelerator driven subcritieal reactor (ADS) works 
well as a solution in that it destroys heavy isotopes 
contained in the used fuel from a conventional nu- 
clear reactor, while at the same time producing elec- 
tricity. The controllable neutron source generated 
by an accelerator also grants the system inherent 
safety. 

In this work, an elementary core layout is cs- 
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Figure 1 : Time evolution of radiotoxic inventory in high-level 
waste (normalized by 1 ton of enriched uranium) 
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tablished for the to-be-constructcd China ADS (C- 
ADS). Based on the core, runtime neutron behav- 
ior, thermal power distribution and transmutation 
performance after long operation time are worked 
out. MCNP5 is applied to simulate neutron behav- 
ior. And a self-developed BUMA code is applied 
to carry out burn-up calculations. The results are 
required for future irradiation damage investigation 
and thermal hydraulics and safety analysis. 

2. Reference core and fuel description 

2.1. Core structure, neutron source and structural 
materials composition 

C-ADS consist of a high-energy proton accelera- 
tor, a Lead-Bismuth Eutectic spallation target and 
a subcritical core. The reactor core, a cylindrical as- 
sembly system, consists of 168 hexagonal prismatic 
fuel assemblies (FA), each of which contains 61 fuel 
elements (FE). Liquid lead-bismuth eutectic alloy 
is used as coolant for its low melting point, high 
boiling point and low neutron absorption cross sec- 
tion. Fuel assemblies and coolant in the core are 
radially and axially covered by a layer of lead re- 
flector. The cross-sectional views of the core, a fuel 
assembly and a fuel element are given in Fig2. 

How a neutron generated from the spallation 
reaction behaves and its dependence on incident 
proton parameters is beyond current investigation. 
Isotropically-emitted neutrons generated from the 
target lateral surface along with those produced 
from fission function as neutron source in the 
model. Energy distribution of neutrons from the 
spallation target is approximated by a Maxwell 
energy spectrum with the temperature parameter 
lOMeV. The FE cladding is made of ferriticmarten- 
sitic steel (FMS) T91 [3], whose typical simplified 
composition is presented in Table 1. 

Table 1: Composition of T91 steel 



Element 


wt% 


Cr 


9.0 


Mo 


1.0 


V 


0.2 


Nb 


0.1 


Mn 


0.5 


Si 


0.4 


Ni 


0.4 


C 


0.1 


Fe 


88.3 




Lead 

Reflector 



(a) Horizontal cross-section of the core layout 




(b) Horizontal cross-section of a FA 

Figure 2: Radial schematics of the 

Core Diameter: 2.21m, Height 2.20m 
Reflector thickness: 0.5m 
FA diagonal: 84.9mm 
FE _Ri:6.00mm, Ro:6.96mm, Hcight:90mm 
Duty Ratio: 49.4% 
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Fuel pins loaded in a FE are assumed to be 
CERMET pellets. They are composite materi- 
als consisting of {Puy,MAi-y)02-x particles dis- 
persed in a metallic molybdenum matrix (called 
thereafter CERMET) [4]. Though in this calcu- 
lation model, plutonium, minor actinides, oxygen 
and molybdenum are homogeneously distributed in 
fuel pins, their atomic fraction and fuel pin density 
are determined from the metal oxides and molyb- 
denum densities and their weight fraction. Sim- 
ulation of dispersed particles will be investigated 
in further research concerning fuel element and 
[Puy.M Ai-y)02-x particle tests. 

2.2. Controlled fuel composition parameters 

Fuel composition optimization is based on the 
isotopic compositions of LWR spent fuels in [5] [6] 
[7] since it is assumed that contemporary technol- 
ogy cycles spent fuels without distinguishing iso- 
topes of an element. For the fixed isotopic com- 
position (given in Table 2), weight fractions of the 
metal dioxides including Pu02, Am02, and Cm02 
are adjusted for better neutron performance. For 
a fuel pellet in the loaded CERMET fuel elements, 
metal dioxide microspheres and the molybdenum 
matrix account for nearly 40 vol% and 60 vol% re- 
spectively. Weight fractions among the metal ox- 
ides and the matrix are optimized for required neu- 
tron performance including subcritical multiplica- 
tion factor ks, effective neutron multiplication fac- 
tor fce//, fiat power distribution, high MA transmu- 
tation efficiency and low fuel irradiation damage. 

3. Computation tools 

The calculations in this work are performed ap- 
plying the probabilistic neutronic code MCNP5 and 
a coupled code system BUMA for rough burning- 
up analysis. MCNP5 yields instant results includ- 
ing effective neutron multiplication factor, neutron 



fiux and various reaction rates averaged over spe- 
cific number of transportation histories. The de- 
veloped burning-up code with nuclear data com- 
piled from evaluated libraries [8] provides an ele- 
mentary view over time-dependent transmutation 
performance. In the calculation progress, MCNP5 
gives spatial reaction rate distribution and BUMA 
uses MCNP5 output file to evolve fuel compositions 
in MCNP5 input files. A batch file calls the above 
two steps alternatively till the end of the transmu- 
tation cycle. 

4. Instant transmutation efficiency and 
power distribution 

The life time of a neutron history is found to be 
not greater than the scale of 10^ shakes. And the 
current, flux, collision reactions accumulated over 
the history can be esteemed as instant values if di- 
vided by the time period or timed by the source gen- 
erating rate. System performance under the above 
equilibrium conditions is believed to be a good ba- 
sis for optimization analysis. And therefore in this 
work, instant transmutation efficiency and power 
distribution are analyzed first and then burning-up 
calculation is performed on the base of equilibrium 
cycles. 

J^.l. Determination of initial fuel composition 

Initial fuel composition aims to result in sufficient 
future MA transmutation as well as maintain safe 
effective multiplication factor k^f / , subcritical mul- 
tiplication factor kg and external neutron source ef- 
ficiency (f*. ks, an important parameter for a sub- 
critical system is defined as below [9]. From the 
neutron balance equation describing a subcritical 
system with an external source (Eq(l)), total fis- 
sion neutrons F (Eq(2)), total source neutrons S 
(Eq(3)) can be determined, and kg defined with F 
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Table 2: Initial isotopic composition 



Isotopic Vector Atom fraction (at%) 

Pu/Np/Am/Cm/Mo/O 4.34/1.84/6.47/1.07/58.91/27.45 

234u/ 235u/ 236u/ 238u 0.0132/0.5793/0.3560/99.0516 
238Pu/ 239pu/ 240p^j/ 241p^/ 242p^j 1.28/62.0/23.5/8.90/4.36 

237Np 100 

Am/ Am/ 243 Am 84.61/0.246/15.14 

243Cm/ 244cm/ 245Cm 2.01/90.01/7.97 



and S (Eq(4)). And the relationship between ks 
and kefT is obtained as: 

[9] . Amount of Am and Cm is raised for more incin- 
eration and that of Pu is decreased to maintain k^j f 
between 0.95 and 1 which helps reduce accelerator 
requirements. As mentioned above, weight fraction 
of the metal oxides and the molybdenum matrix in 
a fuel element pin is controlled as independent vari- 
ables. And the initial fuel composition program is 
presented in Table 3. Volume fraction is estimated 
from the theoretical densities of the metal dioxides 
and molybdenum matrix. And taking into account 
the isotopic composition assumed in Sec. 2. 2, the 
fuel compositions in term of individual nuclides are 
obtained for MCNP input file. 

Tabic 3; Initial materials composition of fuel element pins 



Materials 


wt% 


vol% 


UO2 


13.94 


13.49 


Np02 


1.39 


1.33 


PUO2 


14 


12.92 


Am02 


3.35 


3.05 


CmOa 


2.24 


2.03 


Mo 


65.08 


67.19 



Initial fee/ / for the core loaded above is estimated 
to be 0.963±0.001. kg, the multiplication factor 
with external neutron source is 0.992. And from Eq 
(5), external neutron efficiency (f* is 4.7. In fact, 
since the amount of Am and Cm accounts for less 
than 20% of the total amount of Pu and U, Am 
and Cm accumulation cannot compensate for the 
reactivity decrease resulting from U and Pu deple- 
tion. Therefore this initial loading program main- 
tains stable kef / over operation time and is suitable 
for energy production. 



In comparison, if the amount of Am and Cm ac- 
counts for too much of the fuel, neutron multipli- 
cation factor increases over time due to breeding 
of fissile nuclides like 242 Am. Burn-up analysis sug- 
gests that keff of the system approaches 1.00 within 
a month and such a program is not suitable for en- 
ergy production. On the other hand, for an initial 
isotopic composition in which Am and Cm accounts 
for less than 10% of the amount of Pu and U, the 
duty ratio (fuel volume fraction) of the FAs can- 
not exceed 15% and reactor core efficiency will be 
decreased greatly. 

4-. 2. Runtime neutron energy and spatial distribu- 
tions 

Incineration performance in the core mainly de- 
pends on fiux level and neutron spectrum. In or- 
der to investigate the neutron fiuxes and currents 
as time derivatives, the MCNP-yielding neutron in- 
fiuence and current, which are integrated over the 
cutoff time period and are averaged over all source 
histories, arc timed by the product of the neutron 
yield of the spallation source per incident proton, 
Ygp, and the proton beam current Ip. They are as- 
sumed to me 27.4 n/p and 8.1mA respectively and 
the muhiplier is 1.385 x lO^^n/s. With the multi- 
plier, the results given by MCNP such as tally per 
particle are converted to current and reaction rates. 

Neutron current over the target lateral surface, 
the refiector inner and external surfaces, the sup- 
posed surfaces surrounding 7 zones (The 168 FAs 
in Fig 2 are assumed to be arranged in 7 circles.) 
are given in Fig3. The neutron spectrum shifts to- 
ward high energy from the source surface to the 
inner reflector surface. Moderation of the lead- 
bismuth coolant and energy distribution of fission 
neutrons are responsible to the movement. In com- 
parison with the spectra of source neutrons and 
those through the inner reflector surface, spectrum 
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of the neutron moving out shifts down apparently 
and it is caused by neutron leakage. 

The moderation effect is exhibited more directly 
in Fig3(b), spectrum shifts gradually from source 
to zone I, IIjVII and finally contributes to the ap- 
parent shift of spectrum of neutrons entering the 
reflector or leaving the active zones. 

4-.3. Energy deposition and thermal power distribu- 
tion 

Power distribution can be determined from the 
fuel element heat-depositing reaction rates. Taking 
into account the core symmetry, radial power dis- 
tribution is characterized with the line power den- 
sity along one radius. As given in Fig2(a), the 7 
fuel assemblies located in the same radius are se- 
lected with 5 fuel elements along the same line to 
determine the line power density, which is presented 
in Fig4. The radial power peaking factor is deter- 
mined to be 2.18. Summing up the power density 
of the unique fuel elements, multiplied by repeating 
numbers and FE volume yields the total power, 451 
MW with the assumed Ygp and Ip. 

Besides the radial power distribution and peaking 
power factor, axial heat-depositing rate is required 
for further thermal hydraulics and safety analysis. 
One of the 7 fuel elements in one assembly in the 
4th zone is taken as an example to show the ax- 
ial power distribution. The thermal power profile 
in Fig5 is based on a fuel element divided into 30 
equal- volume horizontal sectors. 

4.4- Minor Actinides transmutation validation 

Energy-dependent neutron current or flux 
through separating surfaces illustrated in Sec4.2 de- 
scribes neutrons traversing the fuel elements and 
structural materials. And their spectra are respon- 
sible for reaction rates spatial distribution. As the 
detected surface moves outward the core, the spec- 
tra shift to low energy while at the same time, re- 
action rates of most concerned types decrease out- 
ward (Fig7). The opposite shift direction of neu- 
tron spectrum and reaction rate can be accounted 
to the fact that neutrons at high energy contribute 
mainly to the concerned types of reaction, namely 
the transmutation of MA elements. 

In addition, microscope fission, absorption and 
inelastic cross-section of MA elements decrease over 
the minimum value and reach a high plain above 
about 3MeV (Fig6). The critical energy is consis- 
tent with the turning point in Fig3. In conclusion. 




10-7 10-G 10-5 10-4 0.001 0.01 0.1 1. 10. 
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(a) Current through spallation target surface, inner and outer 
surfaces of core reflector 
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(b) Current through fuel assembly group Zones I, II, III to VII 
Figure 3: Neutron energy distribution 
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Figure 4: Radial line power density for the core 
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Figure 5: Axial power profile in a fuel element 



though energy spectra shift to low energy in the 
core from center to reflector, reactions induced by 
high energy neutrons account for more weight and 
thus are responsible to the radial reaction decrease. 



neutron fission cross section 




. la-io 10-9 lO-ft ia-7 10-6 lO-? 10-4 0.001 0.01 
energy <mev) 

(a) CTj of Am, Cm isotopes 
neutron fission cross section 




energy (mev) 



(b) <Jf of Pu isotopes 
Figure 6: Fission cross-section of MA elements 

MA transmutation can be further validated with 
the neutron spectrum (Fig8) for the same core lay- 
out but with a thermal neutron source. Radial re- 
action rate decrease is also observed in the refer- 
ence core. However, the energy profiles of neutron 
current through the target lateral surface, zone sep- 
arating surfaces and reflector surfaces differ greatly 
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Figure 7: Transmutation rates of MA elements 



from those presented in Fig3. The spectra shift to- 
ward high energy especially after a similar turning 
point at about leV. Therefore, it can be concluded 
that reactions resulting from neutron at low energy 
predominate in the thermal core and MA transmu- 



tation performance is not promising. 
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Figure 8: Neutron energy distribution of a thermal core 

5. Time-dependent and Burn-up analysis 

5.1. Calculation approach 

Instant transmutation results present nothing 
about liow tlie fuel composition and transmutation 
rate vary with time. Burn-up analysis is necessary 
to confirm whether the obtained results are con- 
served over time and to optimize fuel composition 
and realize sufficient incineration of MA nuclides. 

Evaluated nuclear data file libraries and MCNP 
are applied to obtain transmutation relationships 
(Fig9). Burn- up analysis in this work concentrates 
on MA, their mother and daughter isotopes in neu- 
tron reaction or decay, and fission products that 
have high neutron absorption cross section. Main 
MA isotopes remaining in spent fuel and need in- 
cinerating are "^^"^Np, "^"^^Pu, ^^^Pu, "^'^^Am and 
^^^Cm. With neutron reaction cross section and de- 
cay half-life from IAEA nuclear data service [8] and 
reaction rates from MCNP pre-calculation, other 
nuclides like ^^^Np and ^^^Pu are added to the 
transmutation relationship. Fission products that 
are out of Fig9 but responsible for neutron absorp- 
tion are "^/, "s^e and "SATrf. The resulted neu- 
tron poisons are ^^^Xe and ^'^^Sm. And radioactive 
transformation of the above light nuclides is also 
considered in calculation. 

The temporal change of the atomic densities of a 
related nuclide dNj for a discrete time interval At 
is evaluated in Eq(6), where Nm,i indicates atomic 

density of nuclide m during time interval i, cr^^ in- 
dicates micro cross section breeding nuclide j from 



reaction b of nuclide m, a - J indicates micro cross 
section depleting nuclide j from production of nu- 
clide n in reaction d, a|"2 indicates type b decay 
rate of radioactive transformation from nuclide I to 
j, and A^'^j is type d decay rate of radioactive trans- 
formation from nuclide j to k. 
Take ^'^'^Cm for example, 

■^c'miL^Cm.4,i. Traditional burn-up calculation 
works out cross-section as group constants first and 
then the constants are applied for reaction rates 
determination. In this research, reaction rates of 
several concerned types are obtained from MCNP 
output files and group constant calculation is lim- 
ited to fission yields. Therefore, calculation amount 
may be increased since Monte- Carlo reaction counts 
calculation is substituted for group constants deter- 
mination but burn-up analysis precision can be ob- 
tained from Monte-Carlo results. Burning-up anal- 
ysis is executed by alternative calling MCNP5 and 
the coupled codes system BUMA. Initial fuel iso- 
topic composition and reaction rates along with the 
transmutation flow (Fig9) suffice to determine the 
isotope composition after the burn-up time step. 
MCNP5 works out reaction rates averaged over all 
the neutron histories of concerned types described 
above for the initial core composition. BUMA con- 
verts the MCNP5 results to time derivatives with 
special Ysp and Ip values and performs fuel compo- 
sition variation (Eq(6)). Continuing the above two 
steps will yield fuel compositions, neutron energy 
profile through concerned surfaces and cells, spatial 
thermal power distribution and MA transmutation 
rates as time-dependent variables. 

5.2. Burning-up results 

For the core with initial fuel composition as given 
in Sec4.1, the above burning-up analysis is per- 
formed for 365 days of stable operation. The tem- 
poral evolution of atom density for concerned MA 
isotopes is presented in FiglO. And based on the 
evolved EE composition, kgff and ks is calculated 
to confirm reactor core safety. Over the 365d op- 
eration period, k^ff decreases gradually from 0.963 
to 0.939 and kg remains approximately 0.994. 

Two fuel elements from Zone II and Zone 
VII respectively are selected for evolution data. 
Fig. 10(a)10(b), 10(c) describe the evolution in the 
two EE of HLW isotopes given in Figl. As ex- 
pected from decreased reaction rates toward periph- 
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ni.b 



n.d b,l d,k 



eral FAs, outer FE experiences deeper MA trans- 
mutation. On average, a fuel assembly in Zone II 
will transmutate 245.2g of ^^i^m, 35.9g of '^^^Am, 
110. 8g of '^^'^Np in a year. And at the same time one 
in Zone VII will incinerate 64. 5g of '^'^^Am, 10. 3g of 
'^'^^Am, 29. Ig of jVp. The burned MA amounts 
depend on initial fuel composition. In addition, in- 
evitable accumulation of other nuclides like ^^^Am 
and 2'*^Cm is given in FiglO(d). 

5.3. Refueling of the core after incineration 

Due to the decreasing radial reaction rate profile, 
fuel assemblies in peripheral zones of the core fail 
to experience sufficient transmutation and do not 
deserve subsequent procedure to separate nuclides. 
Therefore, the peripheral assemblies are moved in- 
ward and fresh assemblies are refueled from outside 
to realize deep incineration. After 365 days of evo- 
lution, 36 assemblies in the first 3 zones are replaced 
with 36 ones from the 7th zone. The other assem- 
blies from the last 4 zones can be installed in a new 



core. Besides, the refueling program provides an- 
other approach for flattening thermal power profile 
in addition to applying more than one spallation 
source in the system [10], adjusting inert matrix 
ratio [11] and loading different types of fuel for dif- 
ferent reactor core zones [12]. Thermal power for 
the refueled reactor is 415MW if proton beam re- 
mains constant but the power peaking factor fails 
to decrease obviously because the fuel difference 
caused by annual transmutation cannot make up 
for the spatial neutron spectrum difference. Effec- 
tive neutron multiplication factor keff for the re- 
fueled system is determined to be 0.96, which is 
approximately equal to that of the initial core and 
the reactor is equivalently restarted for a new cy- 
cle. Without the refueling program, keff decrease 
can also be compensated by adding burnable poi- 
sons when starting the reactor [13] and burnable 
poisons utility will be confirmed for C-ADS after 
initial fuel composition and refueling programs are 
determined. 
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Figure 10: Fuel composition evolution 
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6. Summary and conclusion 

A neutronic analysis of the simulated C-ADS re- 
actor core with the composite CERMET fuel with 
Mo matrix, selected for MA transmutation, was 
performed aimed at optimizing the transmutation 
and neutronic performances. Runtime neutron en- 
ergy profile and spatial distribution is presented. 
Thermal power of the core is determined to be 
451MW with the radial power peaking factor 2.18. 
Neutron spectrum shift along with reaction rate 
variation validates the transmutation efficiency of 
MA elements. With an approximate burning-up 
analysis method, fuel composition evolution is pre- 
sented for the 365d operation period. 
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